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Abstract 

We revisit the phenomenology of the Constrained MSSM with right-handed neu- 
trinos (CMSSMRN). A supersymmetric seesaw mechanism, generating neutrino 
masses and sizable lepton flavour violating (LFV) entries is assumed to be opera- 
tive. In this scheme, we study the complementarity between the 'observable ranges' 
of various paths leading to the possible discovery of low energy SUSY: the reach of 
the Cern Large Hadron Collider (LHC), the quest for neutralino dark matter signals 
and indirect searches through LFV processes. Within the regions of the CMSSMRN 
parameter space compatible with all cosmo-phenomenological requirements, those 
which are expected to be probed at the LHC will be typically also accessible to 
upcoming LFV experiments. Moreover, parameter space portions featuring a heavy 
SUSY particle spectrum could be well beyond LHC reach while leaving LFV searches 
as the only key to get a glimpse on SUSY. 



1 Introduction 



This is going to be the decade where we should be able to establish whether low en- 
ergy supersymmetry (SUSY) exists or not. We have three main roads to gain access to 
SUSY: (a) Direct SUSY searches at hadron colliders (b) Indirect SUSY searches in rare 
FCNC (Flavour Changing Neutral Currents) and/or CP violating processes (c) Direct 
and indirect SUSY Dark Matter (DM) searches. Given that low energy supersymmetry 
is rather a framework than a predictive model, one has to assume an explicit realisation 
of SUSY breaking in order to make predictions concerning the above mentioned three 
roads to SUSY. The Constrained Minimal SUSY Standard Model (CMSSM) represents 
the 'prototype' of low energy supergravity having the minimal number of free parame- 
ters and being successful in not departing too violently from the Standard Model (SM) 
predictions. In view of the extraordinary agreement of the SM with all the experimental 
data on the above (a) and (b) points, it is of interest to consider models with features 
approaching those of the CMSSM. 

The CMSSM, in its original formulation, is however not complete, and it has to be sup- 
plemented with some mechanism to generate non-zero neutrino masses and mixings. The 
super symmetric seesaw mechanism represents a very attractive candidate to accomplish 
this task. Hence we are led to study a CMSSM with right handed neutrinos (CMSSMRN). 
We consider this as the 'minimal' low-energy supergravity-inspired setting with massive 
neutrinos. Compared to the CMSSM, the main novelty in the CMSSMRN is the presence 
of lepton flavour violation (LEV) [T^. Information from neutrino masses is nonetheless 
not sufficient to determine all the seesaw parameters [2j, which are crucial to compute 
the relevant LFV rates. To remedy this, either a top-down approach with specific SUSY- 
GUT models and/or fiavour symmetries IH El 1^ or a bottom-up approach with specific 
parameterisations of low energy unknowns have been adopted in the literature [Zl IHl E] ■ In 
the present work, we follow a hybrid approach making the minimal assumptions required 
to successfully predict lepton fiavour violation. Assuming one of the neutrino Yukawa 
couplings to be as large as the top Yukawa, and varying the other relevant seesaw pa- 
rameters, we move from situations where LFV rates are significantly large, and already 
constrain the CMSSMRN parameter space (what we call the ^^best-case" scenario) to cases 
where a significant experimental improvement would be needed to detect fiavour violation 
in the lepton sector {^^ worst- case" scenario). 

In the light of this setting, we can now revisit the above three roads leading to SUSY 
discovery. Namely, we wish to conduct a quantitative analysis taking into consideration 
(a) direct search limits for SUSY at LEP and hadron colliders (b) indirect searches in rare 
FCNC processes, particularly concentrating on LFV, and finally (c) the neutralino relic 
density, which should respect the upper bound on the cold dark matter (CDM) content 
of the universe. We outline a strong complementarity between the three mentioned roads 
for SUSY discovery. We generically find, in the "best case" scenario, that the inclusion of 
the LFV constraints significantly restricts the already "Constrained" MSSM parameter 
space. In particular, in the most favorable cases, if stau coannihilations are responsible for 
the suppression of the neutralino relic density, SUSY might 'appear' in LFV experiments 
even before the advent of the LHC, which will thoroughly probe these regions. In the 
other two parameter space areas allowed by the stringent WMAP upper bound [101 
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the CDM content, namely the A-pole funnel and focus point regions, LFV sensitivity is 
found to be usually larger than the LHC reach, provided the unknown leptonic mixing 
angle Ues is not too small. Indeed, given that, particularly in the focus point region, 
SUSY masses are quite heavy, it might well be that LFV constitutes the only tool at 
our disposal to get SUSY signals. Therefore, in the "best-case" situation, the inclusion 
of massive neutrinos has critical implications on the discovery paths for SUSY. On the 
other hand, in the "worst-case" scenario, LFV experimental studies will at best only be 
complementary to LHC direct SUSY searches. In this case, non-zero neutrino masses do 
not have a significant impact on the CMSSMRN phenomenology, which would be almost 
indistinguishable from that of the CMSSM. 

Some interesting analyses of the impact of LFV on the CMSSM parameter space al- 
ready exist: for instance, in Ref. ^T] the correlation between DM and LFV was considered 
in a bottom-up approach. Recently, in Ref. |jl2j the connections between rare processes 
and seesaw generated LFV in CMSSM have been studied. Our analysis simultaneously 
takes into account high energy collider physics and low energy flavour constraints as well 
as the latest information on the DM content of the universe. Moreover, the hybrid ap- 
proach on the seesaw implementation that we adopt here differs from what was done in 
previous works. We detail our results in the following Sections. 

In the next Section, we review the CMSSMRN model and present our assumptions. 
We also briefly review the present status of 'discovery limits' of SUSY within the three 
discovery methodologies. In Sections 3 and 4, we provide full details on our numerical 
results, mostly concentrating on the "best case" scenario. In Section 5, we stress the role 
of Ue3 in connection with direct and indirect SUSY searches. After a brief summary of 
the results, we conclude with an outlook. 

2 The CMSSMRN and LFV 

The CMSSM is based on minimal supergravity with universal soft terms at the GUT scale, 
Mgut ~ 2 ■ 10^^ GeV. The SUSY particles spectrum at the weak scale is determined, 
through Renormalisation Group (RG) evolution, by specifying the following high scale 
parameters: mo, the common scalar mass, Aq, the common trilinear coupling and M1/2, 
the universal gaugino mass. The remaining low energy parameters of the model are 
tan/3, the ratio of the Higgs vacuum expectation values (vev) and the sign of fi. The 
superpotential of the CMSSM incorporating the seesaw mechanism jTSj can be written as 

W = Wy^ + h^.Ue^^H, + /ir,L.z/|i/2 + Mn^^^vt^,, (1) 

where the leptonic part has been detailed, while the quark Yukawa couplings and the 
parameter are contained in Wyq- i,j are generation indices. represents the (heavy) 
Majorana mass matrix for the right-handed neutrinos. Eq.(P3) leads to the standard seesaw 
formula for the (light) neutrino mass matrix 

= -h'^M^^h" ^vl, (2) 

where is the vacuum expectation value (VEV) of the up-type Higgs field. Hi. Under 
suitable conditions on W and M^. the correct mass splittings and mixing angles in A^;^ 
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can be obtained. Detailed analyses deriving these conditions are already present in the 
literature ^1]. Additionally, the decoupled heavy right handed neutrinos also leave their 
imprints on the soft sector in terms of flavour violating entries P which can be probed 
in LFV experiments. 

The amount of lepton flavour violation generated by the SUSY seesaw crucially de- 
pends on the flavour structure of and M/j, the new sources of flavour violation which do 
not appear in the CMSSM. This dependence can be clearly seen in the Renormalisation 
Group-induced entries in the slepton mass matrices at the weak scale: 



where Mr is the scale of the right handed neutrinos. From above it is obvious that if 
either the neutrino Yukawa couplings or the flavour mixings present in are very tiny, 
the strength of LFV will be significantly reduced. Given that there are several possible 
regions in the seesaw parameter space which may generate the observed mixing in the 
neutrino sector, one has to resort to some assumptions on h'^ and on Mr to make a 
quantitative analysis of LFV within this model. 

In the present work, we take the view that the presence of non-zero neutrino masses 
modifies the CMSSM in such a way that the operative seesaw mechanism significantly 
maximises LFV leading, in the near future, to another sound discovery road for SUSY. 
To this extent, let us make the following assumptions [T^ lUl ITBj: 

(a) At least one of the neutrino Yukawa couplings is of 0{1). For definiteness, we choose 
it to be the third eigenvalue, h^, setting it equal to the top quark Yukawa, ht. Assuming 
neutrino masses are hierarchical, this automatically sets the largest eigenvalue of Mr to 
be rather close to Mqut, around 10^^ -7- 10^^ GeV. The second largest and the smallest 
eigenvalues /i2i can be left unspecified as far as the dominant contribution to lepton 
flavour violation is concerned. However, for a complete neutrino mass model, one needs 
to specify these eigenvalues too. 

(b) The matrix which diagonalises the product h'^h'^ ^ has either a CKM-like structure with 
small mixing or a MNS like structure with large mixing. Notice that the mixing present 
in h'^h'^ ^ determines the amount of LFV, as given by Eq.Q. Together with assumption 
(a), the case with small CKM-like mixing gives a "worst case" scenario, whereas the case 
with large MNS-like mixing provides a "best case" for LFV. 

As far as neutrino masses and mixings are concerned, assumptions (a) and (b) have 
been extensively studied in the literature, and it has been shown that they may lead to 
phenomenologically viable models for neutrino masses jHj. The relevant references can be 
found in the reviews listed in Ref . ^3] . Given these assumptions, the LFV mass-insertions 
at the weak scale in the CMSSMRN are given, for the worst and best case respectively 
as: 




(3) 
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Worst case: 



K)3^ « _M±^A?VW'„ln^ + 0(A;)^ (5) 

^3 

(m^)3i ^ \Vttytd'^Ti^^ + 0[h^) . (6) 



Best case: 



"^l)2i ~ ^2 f^tUesU^sln— + C(/l2), (7) 

^3 



("1)31 « -^^^|±^;»?c/.3c/.3i.i%2: + o(ft3^ (9) 

^3 

where V is the CKM matrix and U is the leptonic 'PMNS' mixing matrix. It is now 
obvious that the seesaw generated mass insertions above lead to various LFV processes, 
hke rare leptonic radiative decays, —>■ e conversion in nuclei J7j, r — 3/i ^Hl- One of 
the features which stands out is that in the best case scenario, the amplitude for any LFV 
process involving /x — e transitions depends on the neutrino mixing angle, Ues, of which 
very little experimental information is known, except for an upper bound (See for ex: 
Ref. jl]). The same statement is also true for r — e transitions. The r — /i transitions 
are instead t/es-independent probes of SUSY, whose importance was first pointed out 
in Ref. JHI- Otherwise, the leptonic flavour violating transitions are now completely 
determined only in terms of the CMSSMRN parameters. 

To appreciate the phenomenological impact of these rare transitions, we list below the 
present and upcoming experimental limits on the Ij — > /j,7 decays. 

Present limits: 



BR{^ 67) < 1.2 X 10-" 1201 

BR{T^fx-f) < 3.1x10-^ ini 

BR{t 67) < 3.7 X 10-^ 1221 

Upcoming limits: 



BR{fi^e-f) < 10-13^10-1^ [23 
BR{T^^i-i) < 10-^ [27 
BR{t 67) < 10-« |22i 

In the next Sections, we will apply these bounds to constrain the CMSSMRN param- 
eter space. 
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2.1 The Three Roads 



The seesaw mechanism operative within our model does not modify the SUSY mass spec- 
trum of the CMSSM significantfy. Thus except for LFV and other rare processes, the 
phenomenology of the CMSSM is very similar to that of the CMSSMRN. Let us briefly 
summarise the status of the mentioned three roads and spell out the various constraints 
we consider in the present work. 

(i) Collider Searches: Direct collider searches will give conclusive evidences for the 
existence of low energy supersymmetry. The searches conducted at LEP already provide 
strong bounds on the CMSSM parameter space. The main constraints are 

(a) The SM-like light Higgs boson mass: 

TUh > 114.1 GeV, and (10) 

(b) The lightest chargino mass: 

m + > 103 GeV. (11) 

These bounds significantly constrain the CMSSM parameter space, especially in the low 
tan (5 and low mass spectrum {^^bulk ") region. Soon, more extended CMSSM regions 
will also be probed by the Fermilab Tevatron^. A by far larger portion of the CMSSM 
parameter space will be probed at the CERN Large Hadron Collider (LHC)^. 

(ii) Rare processes: Low-energy supersymmetry can also manifest itself in various 
flavour and CP violating processes (see [2Zj for a review). In the CMSSM, where for 
most of the processes there is no significant deviation from the SM predictions, the main 
constraint on the parameter space stems from BR{b s, 7) |28j (see also j2H] for recent 
progresses). At present, combining the theoretical and experimental uncertainties, this 
constraint is given by pUj : 

2.16 < BR{b 57) X 10^ < 4.34 (12) 

The other main constraints come from the SUSY contributions to the anomalous magnetic 
moment of the muon, 5a^, and from the Higgs mediated Bs — > /i'''/^"- It is known that 
the anomalous magnetic moment of the muon is capable of putting severe constraints on 
the supersymmetric parameter space. However, at present there are large theoretical and 
experimental uncertainties involved in evaluating the SM contribution to 5a^. The other 
significant constraint might come from the upper limit on the BR{Bs A^^/w ) from 
Tevatron but again as far as the CMSSM is concerned, BR{b 5,7) is mostly 

-'^It has been recently shown that parameter space up to mo < 200 GeV and M1/2 < 250 GeV could 
be probed for a low tan (3 value assuming that the integrated luminosity of the Tevatron reaches 25 fb~^ 

m 

^Here, in terms of mo and M1/2, the maximal region covered by LHC, at an integrated luminosity of 
~ 100 fb~^, will be up to 1400 GeV in M1/2 for mg < 1 TeV. At larger mo > 5 TeV, the maximal reach 
in Ml/2 will instead be around 700 GeV 
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a much stronger bound. We postpone the consideration of these effects to a future work. 

(iii) Neutralino Relic Density and Dark Matter Searches: In the CMSSM 
the hghtest neutrahno Xi is the hghtest supersymmetric particle (LSP) over a very wide 
range of parameters, providing an ideal particle candidate for cold dark matter, xi is 
practically always, to a high degree of purity, a bino, except for a thin strip close to the 
parameter space area where radiative electroweak supersymmetry breaking conditions are 
not fulfilled, named focus point region [33j. 

The current "state-of-the-art" of supersymmetric dark matter studies faces on the 
one side an impressive progress in observational cosmology, which recently led to an 
unprecedented accurate determination of the amount of cold dark matter in the universe 

cni: 

QcDuh' = 0.1126^^09, (13) 

where I^cdm represents the cold matter density of the universe and h parameterises the 
present Hubble rate, Hq = lOO/i Km sec-^Mpc"\ with h = 0.72 ± 0.08. On the other 
hand, extremely accurate numerical packages allow to compute, with a typical accuracy 
of 1% or better, the neutralino thermal relic abundance for a given supersymmetric setup 

We will require that neutralinos either totally or partly contribute to the inferred 
dark matter density. We point out that in most models of low energy SUSY it turns 
out that the computed amount of thermal relic neutralinos is typically larger than what 
indicated in Eq. (fT!^ above^. With this in mind, and allowing for the presence of extra 
DM components beside neutralinos, we will only consider as a constraint the 95% CL. 
upper bound on the thermal relic abundance of neutralinos derived from Eq. ()13|) : 

n^,h^ < 0.129. (14) 

As regards prospects for direct and indirect detection of neutralino dark matter within 
the constrained MSSM, |^ HH] it has been shown that the most favorable parameter 
space points lie in the focus point-hyperbolic branch region. For instance, regarding the 
muon flux from neutralino pair-annihilation in the center of the Earth or of the Sun, one of 
the most promising indirect detection strategies, planned neutrino telescopes |^ US] 
will probe, within the CMSSM, practically only regions at low neutralino masses and 
large higgsino content. On the other hand, it has been shown [H] that direct detection 
spin-dependent searches have typical projected sensitivities lying orders of magnitude far 
from the expected signal. As regards, finally, direct detection through spin-independent 
interactions, once again promising signals are predicted only in the low neutralino masses 
and/or large higgsino fractions parameter space points of the CMSSM. To summarise, at 
present and future facilities, direct and indirect dark matter searches will compete with 

■^It should be noted that in several examples, the presence of extra components in the Universe energy 
density at neutralino decoupling, for instance a scalar field in a Brans-Dicke- Jordan cosmology [HZIj a 
quintessential field |38l I39L BO] , or the shear energy density associated to primordial anisotropics |37| . 
only strongly enhances the predicted neutralino abundance: for instance, the presence of a quintessential 
scalar field may lead to enhancements in the neutralino relic density up to six orders of magnitude |40| . 
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the LHC only in a very limited CMSSM parameter space region, in the focus point strip 
at moderate values of M1/2 , and therefore at low neutralino masses. 

Combining all the constraints given by Eqs.lfTUl ITTl [T^ and in particular Eq. (fT^ . 
apart from the well-known bulk region, at very low soft breaking masses, which will be 
thoroughly probed by Tevatron and the LHC, three main regions survive in the CMSSM 
parameter space: 

• The Stau Coannihilation Regions: If the lightest stau and the lightest neutralino are 
quasi degenerate in mass, efficient stau-stau as well as stau-neutralino (co-)annihi- 
lations [IB] significantly contribute in suppressing the thermal xi ^'^lic abundance^. 
Most of these regions will be accessible to the LHC, at any value of tan /?, with the 
possible exception of very large values of M1/2. 

• Funnel Regions: In these regions, the bino-bino annihilation cross section is greatly 
enhanced through resonant s-channel exchange of the heavy neutral Higgses A and 
H. The conditions required for the resonant enhancement are a large value for tan (3 
and /i < 0, which are necessary to fulfill the relation 2m^o ^ itla- Here LHC might 
cover the parameter space corresponding, roughly, to M1/2 < 1 TeV. 

• Focus point or Hyperbolic Branch Regions: These regions are narrow zones corre- 
sponding to large values of mo, yielding a low value for the fx parameter. The focus 
point region lies close to where radiative electroweak symmetry breaking (EWSB) 
is not valid, and is also relatively unstable numerically. The main feature is that 
a non-negligible higgsino fraction in the lightest neutralino is produced. The LHC 
reach of this region is pretty limited, owing to the large values of rriQ and M1/2, 
yielding very heavy gluinos and squarks. 

The 'smallness' of the surviving regions in the CMSSM is a result of the strict uni- 
versality assumptions on the soft masses at the high scale within the CMSSM. It has 
been shown that, relaxing some assumptions about the Higgs [IHllin] or the sfermion [SHI 
soft SUSY breaking mass universality, other coannihilation partners may arise, such as 
the sneutrino [121120] or the bottom squark |50j, leading to larger allowed regions in the 
parameter space compared to the present situation. 

We compute the neutralino relic density through the darkSUSY package [31], which 
merges the accurate numerical SUSY RGE's evolution package ISASUGRA [HI] with the 
state-of-the-art technique to treat the thermal averaging of cross sections, coannihilations, 
resonances and threshold effects [^. To evaluate the lepton flavour violating rates we 
use the complete formulseof Ref. ^H]; in the mass-eigenstate basis. 

^Similar mechanisms with the Ughtest stop are possible in regions with very large trilinear scalar 
couplings [23 E] • We do not consider these regions in the present work. 
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3 The Canonical (mo - M1/2) Plane and LFV 

We present our results for the effects of LFV on the CMSSMRN parameter space in the 
(mo — M1/2) plane, as customary in CMSSM studies. We then concentrate on the regions 
of the parameter space which satisfy all the cosmo-phenomenological constraints outlined 
in the previous Section, which are plotted as green contours. As a starting point, in this 
Section we present the results for the 'best case' scenario in the plane (Mi/2,'"^o) for a 
fixed value of tan/5 and = 0, and a given sign of /x. We will always show in what 
follows both a small and a large tan /3 case, respectively setting tan /? = 10 and 50 (for 
thourough discussions about the tan/3 dependence in LFV processes see Ref. |21 ITH]). 
while we devote sec. l3.1l to the question of a non- vanishing trilinear coupling Aq. 

In Fig. n we show the isolevel curves of constant BR(r — >■ /i7) for tan /? = 10 and 
> 0. From the figure, we see that the present experimental limit already starts probing 
the region of the parameter space at low mo and M1/2, which is likely to be tested at 
Tevatron. The green patch representing the coannihilation region intersects the isolevels 
corresponding to BR(r — > 7x7) between 10~^° and 10~^^, a level of sensitivity likely to be 
achieved only at dedicated r-factories. In this low tan (3 coannihilation region, which will 
be entirely probed at the LHC, LFV does not significantly impact on the quest for SUSY. 

The situation changes however for the large tan/3 region, where one expects a large 
enhancement in the decay rates (which roughly scale as tan^ (3). In Fig|21 we show BR(r 
fi'y) iso-contours for tan (3 = 50. Notice that we have now chosen /i < so that the A-pole 
funnel regions (see figure caption) also appear within the allowed parameter space. Here 
we see that even the present experimental bound on BR(r — >• fij) can exclude values of 
mo up to 800 GeV and up to 300 GeV in M1/2. The intersection of the isolevel curves with 
the coannihilation region and with the funnel takes place for BR(r — > fi'j) ^ 10^^ and 
extends beyond 10~^, well within the reach of present as well as proposed B-factories. 

The amplitude of —>■ e'j is Ues dependent, as mentioned earlier. In Fig. El we plot 
isolevels corresponding to the present experimental limit on BR(/i — > ej) for various 
values of Ues- Notice that we use here as a parameter the mixing angle S13 = e~''^Ue3. We 
choose again tan/3 = 50 and /i < 0. The results are remarkable. Taking Si3=0.2, close to 
the present experimental upper bound from CHOOZ, would rule out the whole parameter 
space. This fact can be interpreted either as a negative result for a large Ug^ or for a large 
value of tan/3. In any case, it implies relevant consequences for supersymmetric seesaw 
models. In the plot, we show values starting from S13 = 0.05. We see that a value of S13 
ten times smaller than the present bound (si3 = 0.02) can still exclude M1/2 < 750 GeV 
for any value of mo. 

Finally, we show in Fig. E] the isolevel curves corresponding to various values of 
BR(yU —>■ 67) in the worst case scenario, which, we recall, features a CKM-like mix- 
ing in the neutrino Yukawa couplings. Though tan/3 is fixed to a large value, 50, LFV 
rates are nonetheless rather suppressed: the current experimental bound, for instance, 
does not give any constraint on the SUSY parameter space. On the other hand, even 
in this case an improvement in the experimental sensitivity on BR(/i — 67) may lead 
to large part of the cosmo-phenomenologically viable regions. From this point 

on we will however only consider the best case scenario, because it typically gives more 
interesting constraints and elucidates the possible relevance of Ue3 for LFV. 
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Figure 1: The area shaded in green shows the parameter space, in the (Mi/2,mo) plane, 
allowed by all phenomenological and cosmological constraints, for tan (3 = 10, fi > and 
Aq = 0. The lines correspond to various isolevel curves at BR(r ^ fi-f) =3 ■ 10"^ 6 ■ 10'^, 
10^^, 10~^, 10"^*', lO^"*^^. The starting point of each line is dictated by the condition 
m^-^ = nif, which fixes the lowest possible value of mo. The isolevel curves terminate 
where radiative electroweak SUSY breaking (EWSB) conditions cannot be fulfilled, there- 
fore giving the largest possible value of niQ . 
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Figure 2: The area shaded in green shows the parameter space, in the {Mi/2,rno) plane, 
allowed by all phenomenological and cosmological constraints, for tan (3 = 50, /i < and 
Aq = 0. The upper area at large mo represents the funnel region. The lines correspond to 
various isolevel curves at BR(r — > /i7) = 3 ■ 10"'', 6 ■ 10"^, 10"*^, 10~^. As in Fig. [7J each 
line starts at rriy^^ = rrif and ends where REWSB conditions are no longer fulfilled. 
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Figure 3: The area shaded in green shows, as in Fig. the allowed parameter space, in 
the (Ml/2, mo) plane, for tan/5 = 50, fi < and Aq = 0. The lines correspond this time 
to the exclusion curves dictated by the current experimental bound on BR(/i — > cy), for 
various values of Si^ = 0.05, 0.02 and 0.012.- points lying to the left of these curves are 
henceforth ruled out by the current bound BR(/i —>■ cy) < 1.2 ■ 10~^^. The extreme case 
Si3 = 0.2 would rule out the entire parameter space, and is not shown. Again, as in the 
preceding figures, each line begins at m^-^ = nif and ends where EWSB conditions are no 
longer fulfilled. 
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Figure 4: Different isolevel curves of BR(/i — > c-f) for the worst case scenario, i.e. cor- 
responding to a CKM-like mixing in the neutrino mixing matrix, on the same parameter 
space as in Fig. at tan/3 = 50, /i < and Aq = 0. 
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Figure 5: The dependence o/BR(yU cy) on the trilinear coupling Aq, for the particular 
choice fi > 0, niQ = 600 GeV and M1/2 = 300 GeV. The three lines correspond to the 
ratio of BR(/i — > cj) computed at a given Aq and its value at Aq = 0, for three 
different values o/tan/5 = 10, 30, 50. 

3.1 The Dependence on Aq 

Before proceeding to the next section, we pause to comment on the imphcations of a non- 
zero trihnear couphng on the LFV results. As mentioned above, so far the parameter Aq 
has been set to zero. We argue here that any value 7^ would give a larger branching 
ratio and therefore the case Ao = provides a lower bound. 

A non-zero Aq affects the amplitudes of Ij — > lij decays in three ways: (i) The first 
effect comes from RG scaling. From Eq.Q, we see that the RG-generated off-diagonal 
flavour violating mass-insertions are proportional to A^. This is the dominant effect 
coming from a non-zero Aq, and it is independent on its sign, (ii) Another way in which 
Aq affects the branchings is through modifications to the particle spectrum. In particular, 
the parameter fi is sensitive to the value and to the sign of Aq. However, this produces 
only a small effect in the branchings ratios, because /i in the CMSSM is associated with 
the heavier chargino and the heaviest neutralino states. This effect is responsible of the 
manifest, though small, asymmetry for positive and negative values of Aq (see Fig. IS]), 
(iii) Last, the left-right mixing elements of the slepton mass matrix contain Aq in addition 
to the /itan/? factor. This 'left-right' mixing may play an important role in enhancing the 
decay amplitudes fUl 1^2] ■ However, it turns out that for most of the parameter space, 
the main contribution to this mixing is from the /itan/3 terms rather than from non-zero 
values of Aq. 

As a quantitative reference, in Fig El we show, as a function of Aq/tuq, the ratio between 
BR(/i — > 67) and its value at Aq = 0, for different values of tan j3. In agreement with 
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what stated above, the minimum value of the branching is obtained for Aq = 0, and there 
is only a slight dependence with the sign of Aq. Thus all the predictions presented in 
this and the later section have actually to be regarded as lower bounds on LFV rates. 
Let us finally point out that the overall variation induced by a non-zero value of the 
Aq parameter lies within less than one order of magnitude. Finally, we remark that we 
numerically verified that these conclusions do not depend qualitatively on the particular 
parameter space point under consideration. 

4 Coannihilations, Funnels and Focus Point 

LFV processes constitute a meaningful constraint on SUSY models, as the results for 
the 'best case' scenario presented in the preceding section manifestly show. We will 
demonstrate that sometimes they even do better than direct accelerator SUSY searches. 
Given this situation, we now come back to the question of the complementarity between 
the various roads leading to SUSY discovery, directly comparing the LHC reach with the 
prospects for LFV experiments. We therefore now proceed to detail on the impact of LFV 
rates for each of the three allowed regions of the CMSSMRN, for benchmark tan f3 values. 
Namely, we will concentrate on (a) the coannihilation region, (b) the funnel region and 
finally (c) the focus point regions, as mentioned at the end of Sec. 3. 

4.1 Coannihilation Strips 

The defining condition for this region is m^^ ~ m^^ , where fi denotes the lighter stau. In 
the present section, we stick to the points in the (mg — M1/2) plane which saturate the 
limiting condition m^^ = m^^, maximising the extension of the coannihilation strip. We 
further choose two values for tan/3, 10 and 50, setting, in the first case fi > 0, while in 
the second /i < 0. The scalar trilinear coupling Aq is always set to zero (see Sec EH}. 

In Fig. El we show our predictions for BR(r fi'y) and for BR(yU — 67) in the 
tan f] = 10, fi > case. The yellow region dictates the lower bound on the neutralino 
mass, provided by the LEP constraint on the mass of the lightest CP-even Higgs boson 
mh Plj. On the other hand, the cyan region gives the upper bound, dictated by the 
point where m^.-^ = nif^ and Q^-^h"^ = 0.129, i.e. the maximal neutralino mass in the 
coannihilation strip compatible with dark matter constraints ^U]- We also show the 
current and projected experimental sensitivity for BR(/i — 67). We stress that all the 
parameter space points shown in the plot will be within the expected sensitivity of CERN 
LHC, as the latter extends, for this value of tan /3, up to approximately m^^ ~ 550 GeV 
in the coannihilation strip. Clearly, for such a low value of tan LFV rates are rather 
suppressed, and at present one can just exclude a narrow region at low neutralino masses 
provided S13 is close to its present upper bound. Moreover, BR(r — > yU7) lies at least two 
orders of magnitudes below the planned experimental sensitivity. Interestingly enough, 
in case the experimental sensitivity on BR(/i cy) is lowered down to 10~^^, it will be 
possible to detect, within this scenario, fi ^ cy for S13 as low as 10~^. 

The gray shaded band on the BR(r fi'y) line is obtained by varying the parameters 
within the coannihilation region. Although the range of mo at a given M1/2 (and therefore 
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Figure 6: BR(r /i7) and BR(/i — * 67), for various values of S13 = 0.2, 0.05, 0.02 
and 0.012, along the line, in the (Mi/2,mo) plane for tan/5 = 10, fi > and Aq = 0, 
where m^^^ = ntf^, i.e. in the lowest part of the coannihilation strip. The yellow shaded 
area at low neutralino masses is ruled out by the LEP hound on rrth, while points within 
the cyan shaded region at large neutralino masses give an f^^i^^ value which exceeds the 
current WMAP constraint on the cold dark matter density. The region shaded in gray 
indicates the possible values for BR(r — * /i7) which one can obtain varying mo within 
the parameter space showed in Fig. The shape of these shaded regions is analogous 
for the other lines referring to BR(/i 67) 's. We also show the current and projected 
sensitivities for BR(yU — > 67). We stress that all the points showed in this plot are within 
the expected CERN LHC reach at an integrated luminosity ~ 100 fb~^. 
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neutralino mass) is exceedingly tiny, the shaded area is somewhat large, the reason being 
that the iso-level curves of LFV rates, as shown in Fig. ^ are approximately parallel to 
the coannihilation strip. 

As regards the large tan (3 region, we pick the benchmark value tan /? = 50, and choose 
/i < in Fig. [7| (a), again along the coannihilation strip. This time, since the isolevel 
curves intersect the coannihilation area almost orthogonally, the overall dependence on 
the mo spread is completely negligible, and the gray shaded area is vanishingly small (but 
depicted over the r — > /i7 line). The lower bound on the neutralino mass is dictated by 
the upper bound on the inclusive BR(6 — > 57), which for large tan/5 and negative sign 
of strongly limits the low mass region of the parameter space. On the other hand, due 
to stronger couplings in the relevant (co-) annihilation cross sections, the coannihilation 
strip extends up to rather large neutralino masses. In the present case, the reach of LHC 
should approximately coincide with the upper bound on the neutralino mass shown in 
the plot ^26^ . Noticeably, for such a large value of tan /3, the current experimental upper 
bound on BR(/x 67) happens to put severe constraints on S13: we can for instance 
qualitatively conclude that if tan/? is so large, then S13 must be of the order 10^^ or 
less. The S13 = 0.2 line turns out, for instance, to be completely excluded by the present 
experimental bounds. Concerning future improvements on the experimental sensitivity 
on BR(/i — > 67), we notice that all the lines showed in Fig. [3 (a) will certainly be within 
future reach. On the contrary, the situation for BR(r — > yU7) is not equally favorable, not 
even in this large tan/3 scenario. Yet, BR(r fi'-f) of (9(10~^) would start probing this 
region. 

Lastly, in Fig. [7| (6) we show what would happen switching the sign of fi to positive 
values: the lower limit on m^^^ is now given by the lower bound on BR(6 — > 57) , therefore 
excluding a smaller region, and the coannihilation strip is also slightly enlarged towards 
larger masses. The predictions for LFV rates are nevertheless not much affected, except 
for the parameter space which is overall much wider, thus leaving an appealing window, 
at low masses, where LFV processes are particularly large, and even BR(r fi'-f) may 
lie within planned experimental sensitivities. 

4.2 ^-Pole Funnels 

We now proceed in our analysis moving to the A-Pole Funnel regions. The defining 
condition for this region is 2m^j ~ rriA, and again we show (mg, M1/2) points which 
saturate the limiting case for which equality holds. In Fig. |H1 we consider tan /? = 50 and 
^ < 0. Also in this case, though the parameter space, as shown in Fig. |21is by far larger 
than in the coannihilation strip, the spread in the LFV rates is again remarkably narrow, 
as can be inferred from the gray shaded region surrounding the BR(r fi'y) line. As 
before, the lower limit on m^^ is set by the upper bound on the inclusive branching ratio 
b — > 57. The projected LHC reach only extends up to m^-^ < 500 GeV, thus leaving a 
sizable portion of parameter space outside visibility at the future CERN facility. This 
is the first instance where LFV experiments actually compete with CERN LHC as an 
additional road to supersymmetry. In fact, values of < 0.2 are already ruled out by 
the current experimental bounds on BR(/i — > 67). Should the experimental reach for this 
branching ratio be lowered down to 10~^^, we would be able to detect e'j signals 
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Figure 7: (a).- BR(r ft'y) and BR(/i — > 67), for various values of S13 = 0.2, 0.05, 0.02 
and 0.012, along the line, in the {Mi/2,mo) plane for tan/5 = 50, /i < and Aq = 0, 
where m^^ = mf^, i.e. in the lowest part of the coannihilation strip. The yellow shaded 
area at low neutralino masses is ruled out here by the inclusive BR{h s-y) bound, while 
points within the cyan shaded region at large neutralino masses give Q^-^h"^ exceeding the 
current WMAP upper bound on CDM abundance. The current experimental bound on 
BR(;U 67) is shown by an horizontal s^id violet line. The region shaded in gray 
indicates the possible values for BR(r — > /i7) which one can obtain varying rriQ within 
the coannihilation varameter svace showed in the lower areen striv of Fia. IB Also here. 




Figure 8: BR(r —>■ /X7) and BR(/i — > 67), /or various values of S13 = 0.2, 0.05, 0.02 
and 0.012, along the line, in the (M1/2, mo) plane for tan/5 = 50, yU < and Aq = 0, 
where 2 ■ m^^^ = m^, i.e. in the central part of the funnel region. The yellow shaded 
area at low neutralino masses is ruled out by the inclusive BR{b 57) bound, while 
points within the cyan shaded region at large neutralino masses give an f^^i^^ exceeding 
the current WMAP constraint on neutralino relic density. The region shaded in gray 
indicates the possible values for BR(r — »• /i7) which one can obtain varying within the 
funnel parameter space as showed in the upper large green strip of Fig. We also show 
the current and projected sensitivities to BR(/i 67). The expected CERN LHC reach 
at an integrated luminosity ~ 100 fb""*^ is indicated by the vertical orange dotted line: at 
neutralino masses larger than m<^^ ^ 500 GeV LHC will probably not be able to detect 
supersymmetry in the present parameter space setting. 
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in the whole funnel region, provided that S13 > (9(10~^). Interestingly, the large mass 
region lies beyond the expected LHC reach. As in the previous case, we also point out 
that an experimental sensitivity on BR(r — > fi'j) of O{10~^) would already start probing 
this region. 

4.3 Focus Point Regions 

In the focus point region very large values of mo lower the Higgs mixing parameter fi, thus 
entailing the generation of a non- negligible higgsino component in the lightest neutralino^. 
This, in turn, yields an enhancement in the annihilation cross section with respect to the 
pure bino case, together with coannihilation effects with the next to lightest neutralino 
and, more importantly, with the lightest chargino, owing to the mass matrix structure 
of neutralinos and charginos. The combination of coannihilation effects and of a larger 
annihilation cross section forces the neutralino relic density to drop to very low values, 
which may be nonetheless compatible with the current dark matter abundance (see the 
discussion in Sec. 12. ip . 

The focus point region poses several computational problems, since it lies very close to 
parameter space points where EWSB fails, and moreover because it is rather fine-tuned, 
being very sensitive to the input parameters, especially the top mass, rrit. For this reason, 
this region of parameter space has sometimes not been included in CMSSM parameter 
space analyses |1H1- Nevertheless, we include it in our discussion. Indeed, we consider it 
useful to analyze the situation from the LFV rates point of view because it will be very 
hard to probe most of this region at the LHC. 

As before, we once again resort to a low (10) and a large (50) value of tan (3, and we 
choose to show the extreme part of the focus point region, i.e. that at the largest possible 
mo, and hence where the higgsino content is maximal. Our choice is again motivated by 
two considerations: first, the region is sufficiently narrow so that LFV rates along the 
focus point region would appear in any case as single lines; second, choosing the maximal 
possible higgsino content we extend the parameter space line up to the largest possible 
neutralino masses. 

The lower neutralino mass bound in the focus point region, due to the mass vicinity 
between the lightest neutralino and chargino, is typically dictated by the chargino mass 
bound from direct searches. As anticipated, the LHC reach is rather limited in this 
region: it lies around a neutralino mass of ~ 200 GeV. The large masses characterizing 
the sfermion spectrum naturally suppress LFV processes, but nonetheless, even in the less 
favorable case of tan/3 = 10, LFV can probe SUSY up to neutralino masses around 800 
GeV, corresponding to very large values of the soft breaking masses at the GUT scale, 
namely M1/2 ~ 5.5 TeV and ttiq 17.5 TeV. As regards the large tan (3 case, BR(/i — >• ej) 
may be within future experimental reach for neutralino masses in the multi- Te V range, 
provided S13 is of C(10^^ -^ 10^^). Therefore, the focus point region, as well as, in general, 
the case of a higgsino dominated hghtest neutralino, typically tends to favor BR(yU — > 67) 

'"^It should be noted that this situation takes naturally place in several other soft SUSY breaking 
scenarios, for instance with non- universal gaugino masses |53|. where the lightest neutralino is mainly a 
higgsino; relevant consequences for LFV can be found in Ref. |54| . 
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Figure 9: BR(r fi'y)and BR(/i 67), for various values 0/^13 = 0.2, 0.05, 0.02 and 
0.012, along the extreme focus point region, in the {Mi/2,itlq) plane for ta.n (3 = 50, /i < 
and Aq = 0. The parameter space points we use here are those such that the higgsino 
content of the lightest neutralino is maximal. The cyan shaded region at large neutralino 
masses gives an Q^-^h"^ exceeding the current WMAP constraint on CDM density. We also 
show the current and projected sensitivities to BR(/i — > 67). The CERN LHC reach lies 
at neutralino masses smaller than 200 GeV. All the SUSY parameter space points in this 
plot are therefore outside CERN LHC reach at an integrated luminosity ~ 100 fb"^. 
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with respect to direct accelerator searches in the quest for supersymmetry^ . 

5 The role of Ues 

The importance of the parameter 11^3 in its connection with the neutrino mixing matrix 
and lepton flavour violation in the SUSY seesaw is apparent from the discussions we car- 
ried out in the previous two sections. The precise value of Ues turns out to be of critical 
importance in at least two contexts in the present analysis: 

(i) LFV versus LHC 

Within the coannihilation regions, as we have seen, LFV will play only a "supporting" 
role to the more powerful LHC searches. The first instance where LFV may reveal itself as 
a superior tool arises in the heavy mass A-pole funnel region. However, the prominence of 
LFV searches is crucially dependent on Ues- To make this more precise, we plot in Fig. ^2 
the BR(/i — >■ ej) as a function of S13, showing iso-neutralino mass curves. From the plot it 
is evident that an experimental sensitivity of 10~^^ would allow to detect LFV as long as 
S13 > 5 ■ 10~^. On the other hand, there is a large band, at heavy SUSY particles masses, 
lying beyond LHC reach, which will be fully accessible to LFV experiments. The same 
holds true in the focus point region, where, provided f/gs is not too small, BR(yU — 67) 
will probe SUSY far more effectively than the LHC. 

(ii) /i — i> 67 versus r — /i7 

We have seen that, if S13 is not too small, the constraints coming from /i — 67 are stronger 
than those derived from r — > fij. Hence, there should be a critical value of S13 below 
which T ^ fi'-f becomes more relevant than fi —>■ 67. What is exactly this value? To 
answer this question we show in Fig. BR(/i — * 67) in units of 10~^^ as a function 
of Si3, for the particular point at ttiq = 343, M1/2 = 500, tan/5 = 50 and /i > 0. The 
dashed-dotted line denotes the value of BR(r — > fi'y) in units of 10~^. We observe that 
the intersection between the t ^ fij line - which is almost independent of S13 - and the 
IJ, ^ ej line takes place at a value of S13 close to 10~^. If S13 > 10~^ then /i — 67 is more 
likely to be observed in future experiments than r — >■ n'f. Although these results were 
obtained for a specific value of mg and M1/2, changes in these parameters are expected 
to affect equally both processes so that the critical value S13 ^ 10^'^ does not depend on 
mo or Ml/2 (We have also checked numerically that this is indeed the case). When S13 is 
very small, the contribution to BR(/i — > 67), proportional to the second Yukawa coupling 
dominates over that proportional to the top Yukawa coupling. In same figure, we show 
the predictions for BR(yU —>■ 67) with and without taking into account the effect of the 
second Yukawa coupling, which we set equal to the charm Yukawa jB]. Notice that such 
effect is only relevant for S13 < 7 x 10~^. 

^In a recent paper it has been pointed out that at future e+ Linear Colhders with center of 
mass energy ^/s « 0.5 1 TeV, the accelerator reach in the CMSSM focus point region may be by far 
larger than that at the CERN LHC 
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Figure 10: BR(r fi'y)and BR(/i e-y), for various values of S13 = 0.2, 0.05, 0.02 
and 0.012, along the extreme focus point region, in the {Mi/2,rno) plane for tan/S = 50, 
yU < and Aq = 0. Again, the parameter space points we use here are those such that 
the higgsino content of the lightest neutralino is maximal. We show neutralino masses up 
to 3 TeV, which are still allowed by relic density considerations. The expected sensitivity 
of CERN LHC is showed by the vertical orange dotted line, while the yellow shaded area 
on the left indicates the bound stemming from the chargino mass limit set by LEP direct 
searches. 
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Figure 11: The dependence o/ BR(/i cy) on si^, for tan/3 = 50, /i < and Aq = 
along the parameter space line, in the (Mi/2,'^0) plane corresponding to 2 ■ m^-^ = rriA, 
i.e. in the central part of the funnel region. The upper region, shaded in yellow, is ruled 
out by the 6 — > 57 bound, while the lower region, shaded in cyan, is disallowed by the 
VLy.^h'^ bound on the neutralino relic density. The three lines respectively correspond to the 
lower and the upper neutralino mass limits and to the largest neutralino mass within LHC 
reach. We also report the current experimental upper bound on BR(/i — 67) . 
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Figure 12: The ratio of the predicted BR(r ji^y) {respectively BR(/i — 67)) and the 
approximate projected maximal sensitivity of 10~^ {resp. 10""^^) at a particular parameter 
space point along the coannihilation strip at tan/5 = 50, vanishing Aq and positive /i. The 
green dotted line corresponds to the case where the effect of a second non-zero neutralino 
Yukawa coupling, set equal to the charm quark Yukawa coupling, is taken into account. 
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6 Summary of the Results 



Here we summarise our results on the complementarity of the three search roads for the 
three allowed regions of the CMSSMRN in the 'best case' scenario. 

• Coannihilation Regions: In these regions, which are mostly accessible at LHC, an 
improvement of two orders of magnitude in the branching ratio sensitivity would 
make — > 67 visible for most of the parameter space as long as S13 > 0.02, even for 
the low tan j3 region. For large tan (3, independent of S13, t fi'-f will start probing 
this region provided a sensitivity of O{10~^) is reached. 

• A-pole funnel Regions: In these regions the LHC reach is not complete and LFV 
may be competitive. If S13 > 10~^, the future /i — > 67 experiments will probe most 
of the parameter space regions. As before, r ^ will probe this region once the 
BR sensitivity reaches O{10~^). 

• Focus Point Regions: Since the LHC reach in this region is rather limited due to the 
large mo and M1/2 values, LFV could constitute a privileged road towards SUSY 
discovery. As we pointed out in Sec. 121 DM searches will also have in future partial 
access to this region, leading to a new complementarity between LFV and the quest 
for the cold dark matter constituent of the universe. 

7 Outlook 

So far, the three roads towards SUSY we mentioned are (a) direct accelerator SUSY par- 
ticles searches; (b) indirect searches through rare FCNC/CP violating processes and (c) 
direct and indirect DM searches. In this decade, these three ways are going to remain 
the best tools to get evidence for SUSY with foreseeable impressive improvements on 
their sensitivity reach (in particular, for (a), with the advent of LHC). Our work tries 
to provide a critical assessment on their complementarity, focusing on the maximal reach 
of each of them in various viable parameter space regions. To make such an analysis 
quantitative, one needs a specific low-energy realisation. We chose to work in the con- 
text of the CMSSMRN model, because (i) 'CMSSM' is the prototype of 'safe' flavour 
blind low-energy SUSY extensions to the SM compatible with all the precision SM tests 
and (ii) 'RN' provides an appealing mechanism to generate neutrino masses. This model 
is severely constrained by direct SUSY particles searches, FCNC constraints, and, even 
more strongly, by the information we gained about the amount of CDM from WMAP and 
large scale structure analysis. As a matter of fact, only three regions of the CMSSMRN 
parameter space still survive. In the most favorable ('best') case, LFV is certainly com- 
plementary to LHC searches in all these three regions and, interestingly, it can probe 
portions of the parameter space which will be unaccessible to the LHC. 
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